Entanglement entropies in free fermion gases for arbitrary dimension 
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We study the entanglement entropy of connected bipartitions in free fermion gases of A'^ particles 
in arbitrary dimension d. We show that the von Neumann and Renyi entanglement entropies grow 
asymptotically as TV^'^"^-'/'^ In A^, with a prefactor that is analytically computed using the Widom 
conjecture both for periodic and open boundary conditions. The logarithmic correction to the 
power-law behavior is related to the area-law violation in lattice free fermions. These asymptotic 
large- A'' behaviors are checked against exact numerical calculations for A^-particle systems. 



The study of the entanglement properties of many 
body quantum systems has attracted much attention in 
recent years, mainly in relation with the critical behavior 
displayed at quantum phase transitions [ll-Q. In partic- 
ular, lots of studies have been devoted to quantify the 
highly nontrivial connections between different parts of 
an extended quantum system, by computing von Neu- 
mann (vN) or Renyi entanglement entropies of the re- 
duced density matrix pA = Ti'bP of a subsystem A with 
respect to its complement B. 

While a good understanding of the entanglement has 
been achieved in one-dimensional (ID) systems , where a 
leading logarithmic behavior has been established at con- 
formal invariant quantum critical points 04^1 in higher 
dimensions the scaling behavior of the bipartite entan- 
glement entropy is generally more complicated, with- 
out a definite general scenario for the universal behav- 
ior at a quantum phase transition, despite of an enor- 
mous effort on the problem. Earlier investigations in 
the black hole physics proposed the general validity 
of the so-called area law for d > 1 the entangle- 
ment entropy is asymptotically proportional to the sur- 
face of the area separating the two subsystems A and B. 
The area law has been generally proven for gapped sys- 
tems Q , independently of the statistics of the microscop- 
ical constituents. This result is corroborated by exact 
and numerical calculations for systems with topological 
order |^], at 2D quantum Lifshitz pointsjlfl], and some 
explicit field theoretical computations The situa- 

tion is less clear for critical systems. On the one hand, 
the area law holds in bosonic systems 0, [H, [11] . Quan- 
tum critical points described by a Landau-Ginzburg ac- 
tion satisfy the area law , which is also supported by 
holographic calculations by means of AdS / GET [l^ with 
subleading logarithmic corrections. Even in some disor- 
dered systems the area law is confirmed by analytical and 
numerical calculations and the presence of universal sub- 
leading term has been shown [1^. On the other hand, 
the simplest condensed matter models, i.e. free fermions 
on a lattice, show multi plic ative logarithmic corrections 
to the area law [l^ [T7l424| . i.e. for a large subsystem 
A of linear size £ in an infinite d-dimensional lattice the 
entanglement entropies scale like 



where S'*^"^ (A) are the Renyi entropies defined as 

S(")(A) = -^lnTrpS. 
1 — a 
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Eor a — > 1, one recovers the vN definition S{A) = 
S'^^^A) = — Tr [p^lnp^]. It has been argued that these 
logarithmic corrections should also appear for interacting 
fermions with a finite Eermi surface |17h24| . 

In this paper we investigate the entanglement prop- 
erties of free fermion gases of N particles in arbitrary 
dimensions, focussing on the large- asymptotic behav- 
ior of the entanglement entropies of connected spatial bi- 
partitions. Analogously to lattice models, an interesting 
issue concerns the existence of multiplicative logarith- 
mic corrections to the large- A^ power-law behavior corre- 
sponding to the area law. A related question is whether 
these logarithmic corrections are due to the presence of 
corners in the area separating A and B, or they are also 
present for a smooth surface. 

In Ref. [2^ we developed a systematic framework to 
calculate the entanglement entropies of free fermion gases 
for a finite number of particles A, based on the overlap 
matrix of the lowest A^ states. This method has been 
already applied to ID systems [25l - l27j . but it allows also 
computations in higher dimensions. We consider systems 
of A particles in a finite volume of arbitrary dimension d, 
and focus on the large- A^ scaling behavior of the spatial 
entanglement entropies of connected bipartitions, finding 
a general A^('^~i)/'' In A^ behavior for large A. This multi- 
plicative logarithmic correction to the large- A^ power law 
is related to the logarithmic correction to the area law in 
lattice models. Gombined with simple scaling arguments, 
this result allows us to conclude that multiplicative log- 
arithmic corrections to the area law are general features 
of frce-fermion systems both on the lattice and in the 
continuum. 

The method. The ground-state wave function of a gas 
of A noninteracting spinless fermions is 5'(xi, ...,XAr) = 
det[1lJi{xj)]/^/Nl, where -ipi are the normalized wave 
functions of the one-particle problem with lowest en- 
ergies. Explicit expressions of the one-particle eigen- 
solutions are products of eigenfunctions of correspond- 
ing ID Schrodinger problems, i.e. V'iii,n2,- -,nd(^) = 
riiLi "^"i where the nt label the eigenfunctions along 
the d directions which are solution of — (9^0„(x) = 
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en4>n{x), SO that the energy of the d-dimensional problem 

is En-^.n2,...,nd = J2i=i^ni- Note that, although the spa- 
tial dependence of the one-particle eigenfunctions is es- 
sentially decoupled along the various directions, fermion 
gases in different dimensions present notable differences 
due to the nontrivial filling of the lowest N states to ob- 
tain the ground state of the A^-particle system. 

In the following we consider hypercubic systems of 
linear size L with periodic boundary conditions (PBC) 
and Dirichelet (open) boundary conditions (OBC). 
The corresponding normalized ID eigensolutions are 
= i-i/aeS'T^fc^/i for PBC and 0^^^(a;) = 
(2/L)i/2sinfc7rx/L for OBC, with fc e Z and a; G [0,L]. 

The spatial entanglement entropies of connected bipar- 
titions can be computed using the method developed in 
Refs. [2^, [2^ (to which we refer for details). For this 
purpose, we consider the N x N overlap matrix A with 
elements 

A„m = / d''zV^(z)V'm(z), n, m = 1, A^, (3) 

J A 

where the integration is over the spatial region A, and 
involves the lowest TV energy levels. Then, the Renyi 
entanglement entropies are derived from the formula 

JV 

5(")(A) = ^e„(a„), (4) 

where a„ are the eigenvalues of A, and 

ea(A)--^ln[A" + (l-A)"]. (5) 
1 — a 

The vN entropy is obtained by the limit a — >■ 1. The 
knowledge of the 5^"-' (A) for different a characterizes the 
full spectrum of non-zero eigenvalues of pA [11] , and gives 
more information about the entanglement than the single 
vN entropy. Furthermore, numerical stochastic methods 
such as classical [1^ and quantum Monte Carlo [sO] can 
measure only Renyi entropies with integer a. 

For systems with both PBC and OBC, the form of the 
overlap matrix is analogous to that of the two-point cor- 
relation of free spinless fermions on the lattice. Indeed, 
taking as simpler example the PBC case, the two-point 
function is 

Giat(x,y)^/ e'>^ (''--) (6) 

where r(/i) is the volume limited by the Fermi surface 
9r(/i), and x, y arc the positions of the sites of the lattice 
in unit of the lattice spacing, thus x G Z"*. On the other 
hand, after inserting the one-particle eigensolutions, the 
overlap matrix ([3]) reads 

A„„ = / d''x e^2-(k,„-k„).x 

J A 
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FIG. 1: (Color online) The bipartite vN entanglement entropy 
S for space regions A — [0,ziL] x [0, 22^] • • • x [0, ZdL] in 
2D (top) and 3D (bottom) systems with PBC. We report 
S/N^'^^"^ vs A''. The full lines show the asymptotic behavior 
5 = cln iV -I- 6 where c is the geometry-dependent constant in 
Eqs. (|14|l and (|15|l . while the constants b are fitted from the 
data for the largest values of A''. 

Eqs. ^ and ([7]) look very similar and present a clear 
correspondence: in the lattice two-point function the in- 
teger arguments are related to the sites and the inte- 
gration over the momenta is within the Fermi surface, 
while in the overlap matrix the indices are related to the 
states and the integration is within the spatial region A. 
Thus N plays the role of the volume of the subsystem in 
the lattice model. This correspondence has been already 
exploited to derive rigorous asymptotic large- A'^ results 
for the bipartite entanglement entropy in ID fermion 
gases [25l. [26j. 

Systems with periodic boundary conditions. To begin 
with we consider systems with PBC. We compute the en- 
tanglement entropies for regions A — [0, ziL] x [0, Z2L] x 
... X [0, ZrfL]. If all < 1, ^ is an hypercube with cor- 
ners, but if all Zi = 1 except one (let us say zi < 1), then, 
because of PBC, the subsystem A is d dimensional strip. 
In order to calculate the entanglement entropies, we con- 
struct the overlap matrix A from the one-particle wave 



3 



functions, calculate numerically its eigenvalues and then 
use Eq. ^ for S'^°'\ Practically, the data are obtained 
for increasing values of the energy E, and N is extracted 
by counting the states with Ei < E. The results in d = 2 
and d — 3 are reported in Fig. [1] for some values of Zi and 
for N up to O(103). The plot of S/N^-^^'^ vs N clearly 
confirms the presence of the multiplicative logarithmic 
corrections. 

For the lattice free fermion model, the coefficient of the 
leading large-^ behavior in Eq. JT]) has been evaluate d by 
assuming the validity of the Widom conjecture [2l], [3l| . 
In the thermodynamic limit, the bipartite entanglement 
entropy of the region A can be computed from the lattice 
two-point function ^ restricted to A. In the case of a 
connected space region A of size i, the Renyi entangle- 
ment entropies behave as [3 HH 
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where ^ is the chemical potential, D, is the real-space re- 
gion A rescaled by £ so that its volume equals one. n^: 
and rifc denote the normal vectors on the spatial surface 
dfl and the Fermi surface 9r(/i), respectively. The valid- 
ity of these predictions has been supported by numerical 
results [Ulli. 

Calculations based on the Widom conjecture allow us 
to derive an analogous result for the large- asymp- 
totic behavior of the entanglement entropy in fermion 
gases defined in the continuum, as already suggested by 
the above-mentioned correspondence between the over- 
lap matrix A and the lattice two-point function. For a 
d-dimensional cubic volume of size L with PBC, con- 
taining A'^ fermions, the corresponding Fermi surface is 
a d-dimensional sphere of radius kp, which is related to 
the particle number N by 



Vr = vdkf - N, 



rd/2 



Vd 



r(d/2 + i) 



(10) 



We want to compute the vN entanglement entropy of 
a space region A, which may be given by the volume 
[0,ziL] X [0,Z2L] X ... X [0, ZdL] with Zi < 1. Using the 
Widom conjecture as in Ref. |18| . we obtain the asymp- 
totic large- behavior 

S{A) =cN^'^-^'>/'^lnN + o{N^'^-^'>/'^\nN), (11) 



da 



where da is the surface of the rescaled volume [0, zi] x 
[0,^2] X ... X [OjZci], and is the Fermi surface rescaled 
to have a unit radius. The integral ([T^ gives 
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FIG. 2: (Color online) The bipartite vN entanglement entropy 
of space regions A — [0,ziL] x [0,Z2i] in a 2D system with 
PBC. We plot the quantity S/N^^'^ - clniV versus A^"^ for 
square regions (top) and stripes (bottom). The coefficients c 
are given by Eqs. (jldp and H15|) respectively. 
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d^2, 
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^ 6 V6/ 
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{ziZ2 + Z1Z3 + Z2Z3), d = 3. 



Using the Widom conjecture, one can also derive analo- 
gous formulas for generic shapes of the connected space 

region A. Note that the factor 2 Hi^j '^i J'^^* 

area of the surface separating A and B in units of L = 1. 
This could lead to the erroneous conclusion that the area 
law is satisfied in these systems, but Eq. ([Tl]) is just 
an asymptotic expansion in A^ and the size-dependent 
multiplicative logarithmic correction to the area law is a 
subleading term of the form N^~^l''- times the logarithm 
of the area (see the conclusions for a discussion of this 
issue). 

Note that Eq. is not valid when one Zi is one. 
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FIG. 3: (Color online) We plot N''^/'^ {S^'^^ / A) built using 
the a = 2 Renyi and vN entropies of spatial regions A = 
[0, z\L] X [0, Z2L] for some values of zi, zi, versus N~'^ . 



because the terms at the edge do not contribute to the 
area, i.e. the approach to the asymptotic large- behav- 
ior (jlip is not uniform as a function of zi. In particular, 
when all zi = \ except one (let us say zi), the area sep- 
arating A and B is just 2 (in units of V) and we have 



cs 



r(l) 



6r(^)r(f + i) 



(15) 



independently of z\. In this case, the subsystem A is a 
d-dimensional strip without corners (because of PBC). 
Note that the approach to the asymptotic large- be- 
havior (jlip is not uniform as a function of z\ as well, 
because for zi = or 1 the entanglement entropy triv- 
ially vanishes. 

The Widom conjecture does not provide information 
on the corrections to the leading behavior. As we shall 
see, the analysis of finite- A'^ computations with increasing 
N indicates that they are 0(N^~'^l'^), i.e. 



S(^A) = N^'^-^)!'^ [clnA^ + 0(1)] 



(16) 



with c given by Eq. (|12p . Note that setting d — \ 
into this formula, we recover the rigorous ID results for 
the asymptotic entanglement entropy [IE, [2^. However, 
in ID the Fisher-Hartwig conjecture not only gives the 
asymptotic behavior as above, but also the 0(1) term 
and the precise form of all subleading corrections, both 
for particle systems [H, [2^ and lattice models [H, IsS] . 

The asymptotic behaviors predicted by Eqs. (|lim2p 
are definitely supported by calculations of the bipartite 
vN entanglement entropy, as shown by Fig. [TJ To be 
more quantitative, in Fig. [5] we report for the difference 



N-^I'^S[A) -clnA^ 



(17) 



for 2D systems, with c given by Eq. ((T^)- This differ- 
ence appears limited, although presenting sizable oscilla- 
tions which enlarge with increasing z\ and z^. Its large- A'^ 
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FIG. 4: (Color online) The bipartite vN and Renyi entangle- 
ment entropy of space regions A = [0, z\V\ x [0, z-2.V\ x [0, zaL] 
in a 3D system with PBC. We plot the quantity S/N^^^ - 
Ca InN versus A*'^^ for different values of Zi. 



limit provides information on the subleading behavior of 
N~^/'^S, which appears to be 0(1), supporting Eq. (|16l) . 

The amplitude of the oscillations modulating the dif- 
ference (fT7)) appears to decrease. However it is not clear 
how they decrease and whether the amplitude of the os- 
cillations get suppressed in the large- A' limit. Oscillating 
corrections to the leading behavior have been observed 
also for ID systems, but in the case of PBC only for 
Renyi entropies with a ^ 1 [2^ [s^ . Their origin can be 
traced back to the presence of relevant operators localized 
at the borders of the region A 



34| . However, the ones 



observed in d > 1 are present also for a = 1 and their 
structure (periods, amplitudes etc.) appear very differ- 
ent from the ID ones. For these reasons both the general 
structure and their physical origin remains unclear and 
this issue would deserve further investigations. 

For the Renyi entropies, the application of the Widom 
conjecture leads to 
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-S{A)+o{N'^'^-^'^/'^\nN), (18) 



analogously to lattice models [2H. This is supported by 
the data shown in Fig. [3l where we plot the difference 
Af-i/2(S'(2) _ 35-/4) between the a = 2 Renyi and vN 
entropies of spatial regions A = [0, ziL] x [0, Z2L\ for 
some values of zi , 2:2 . The data appear to approach a 
constant but with large oscillations that make difficult 
to extract the rate of convergence. 

In Fig. |4]we report analogous results in d = 3. Again, 
the subtracted quantity A^~2/35'(") — CalnA^ turns out 
to be limited, indicating that it is 0(1) (in agreement 
with Eq. ([T6| ). with sizable oscillating corrections whose 
qualitative features are again unclear. 

Open boundary conditions. The results for the asymp- 
totic large- A^ behavior of the entanglement entropy can 
be extended to OBC. For simplicity, we consider only 2D 
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FIG. 5: (Color online) The half-space vN and a = 2 Renyi 
entanglement entropies for a system with OBC, up to A'' ~ 
3000. We plot the JV-^/^^W _ y^j^ .^^j^-jj -^^ ggj 



systems. We first study the half-space entanglenient en- 
tropy, i.e. when A corresponds to the region zi < L/2. 
The Widom conjecture gives 



-,(q) 

'hs 



7V1/2 
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127ri/2 



lniV + 0(I) 



(19) 



where the coefficient of the leading term is half of the 
one for PBC, due to the fact that the boundary between 
the half spaces is half of the one of stripes. Fig. [5] show 
results at fixed up to w lO'^, for the half-space en- 
tanglement entropies S^^^ showing agreement with the 
prediction and the presence of sizable oscillating correc- 
tions as for PBC. 

The presence of boundaries breaks translational invari- 
ance. Thus the entanglement entropy depends on the lo- 
cation of the subsystem A. We consider the case when 
A corresponds to the region limited by an hyperplane 
parallel to one axis and at distance x + L/2 from the 
boundary. We denote the corresponding entanglement 



entropy by s'"^\x) (with sl^g 
the quantity 



S'^g\0)). Fig. reports 



S^^\x) ^ S^^\x) - S^g"\0) = Ari/2 [/^")(^) + oil) 

(20) 

where the function /a (a;) is a 0(1) contribution and can- 
not be obtained from Widom conjecture. Fig. [6] shows 
clearly that, despite of the oscillating corrections, the 
data for many different A'^ collapse on the same master 
curve supporting the scaling ansatz. 

We also consider the case of a subsystem A enclosed by 
two parallel hyperplanes at the same distance x < L/2 
from the center, both parallel to one of the cubic axes. 
The corresponding entanglement entropies are denoted 
by S's"\x). Using again the Widom conjecture for the 
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FIG. 6: (Color online) Bipartite vN entanglement entropies 
Sa(x) (top) and Ss{x) (bottom) for OBC as defined in Eqs. 
(Uni) and ((2T1l. 



leading large- A^ behavior, we expect 
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67ri/2 



■\nN + /s"\x)+o{l) 
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where the coefficient of the leading term equals the one 
for PBC. Again the functions fs{x) cannot be obtained 
from Widom conjecture, but the validity of this ansatz is 
supported by the data reported in Fig. |6l 

Conclusions. We have studied the vN and Renyi entan- 
glement entropies of a system of A'' particles in a volume 
L'* with periodic and open boundary conditions. The 
large A^ asymptotics of these entanglement entropies can 
be obtained by means of Widom conjecture Jilj and can 
be written in the simple form 



1 



(22) 



where the coefficient cs is given in Eq. p3|) while Az is 
the area between A and B in units where L — I. 

This large- A' asymptotic behaviors can be also turned 
into an asymptotic size dependence of the entanglement 
entropies for connected bipartitions of fermion gases in 
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a cubic box of volume L , in the thermodynamic limit 
when TV, L — > cxj and N / L'^ ~ p with p the particle den- 
sity. Eq. (|22|) can be written as 

-^^^A..csp'-'/'L^-'\nL^-^p'-'/'^A\nA, 
1 + a ^ a— 1 

where we used L^^^Az — A with A is the area sepa- 
rating A and B. This confirms the presence of general 
logarithmic corrections to the area law in free fermion 



gases. 

While the leading asymptotic behavior in N of S'^") is 
understood from our analysis, the numerical data show 
the presence of oscillating subleading corrections to the 
scaling whose structure and origin is still not clear. 
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